gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十五章/eg15_6_2.m

    
%频域法整定
num=1;
den=conv([1 1 0],[1 5]);
G0=tf(num,den);
margin(G0)
[Kc,pm,wcg,wcp]=margin(G0);
[Kc,pm,wcg,wcp]
%%%%%%%%%%%%%%%%%%%%%%%%%

t=0:0.01:25;
num=1;
den=conv([1 1 0],[1 5]);
G0=tf(num,den);
[Kc,pm,wcg,wcp]=margin(G0)
Tc=2*pi/wcg;
Kp1=0.4*Kc;
Ti1=0.8*Tc;
s=tf('s');
Gc1=Kp1*(1+1/Ti1/s);
Kp2=0.6*Kc;
Ti2=0.5*Tc;
Td2=0.12*Tc;
Gc2=Kp1*(1+1/Ti1/s+Td2*s);
step(feedback(G0*Gc1,1),':',t);
hold on;
step(feedback(G0*Gc2,1),t)
hold off;